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O , Abstract. Dynamical interactions that occur between objects in dense stellar systems are partic- 

1/'^ . ularly important for the question of formation of X-ray binaries. We present results of numerical 

' simulations of 70 globular clusters with different dynamical properties and a total stellar mass 

^~r^ , of 2 X 10 Mq. We find that in order to retain enough neutron stars to match observations we 

^ ^ ■ must assume that NSs can be formed via electron-capture supernovae. Our simulations explain 

^ ' the observed dependence of the number of LMXBs on "collision number" as well as the large 

r^ , scatter observed between different globular clusters. For millisecond pulsars, we obtain good 

^Lc agreement between our models and the numbers and characteristics of observed pulsars in the 

JL ' clusters Terzan 5 and 47 Tuc. 
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1. Introduction 



In globular clusters (GCs), neutron stars (NSs) are seen as low-mass X-ray bmaries 
(LMXBs), bright or quiescent, and as binary or single millisecond pulsars (MSPs). As the 
numbers per unit mass of LMXBs and MSPs in clusters greatly exceed their numbers in 
the Galaxy, their origin has been linked to stellar encounters that should occur frequently 
r — ' in an environment of high stellar density (Clark 1975). Galactic GCs are known to contain 
f^ , thirteen bright X-ray sources, and as many as ^100-200 quiescent LMXBs (qLMXBs) 
are thought to exist in the Galactic GC system (Heinke et al. 2003). Millisecond pulsars, 
likely descendants of LMXBs (for a review see, e.g., Bhattacharya and van den Heuvel 
^ ' 1991), also are present in GCs in great numbers: about 140 GC millisecond radio pulsars 
have been detected by no\4il, with more than a dozen in several GCs - in 47 Tuc (Camilo 
et al. 2000), M28 (Stairs et al. 2006), and Terzan 5 (Ransom et al. 2005). As only a few 
X-ray binaries or few dozen MSPs are present per fairly massive (more than a million 
stars) and dense cluster, computationally this problem is very challenging. The target 
time for direct A^-body methods to address the million-body problem is around 2020 (Hut 
2006). In our studies, we use a modified encounter rate technique method, described in 
detail in Ivanova et al. (2005), and with the updates described in Ivanova et al. (2007). 



t See http://www.naic. edu~pfreire/GCpsr. html for an updated list 
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Figure 1. The retention fractions as a function of escape velocity (for stellar evolution unaf- 
fected by dynamics) for a Hobbs et al. (2005) kick distribution. Dotted and dash-dotted lines 
show the retention fractions for single and binary populations, core-collapse NSs only. Solid 
and short-dashed lines show the total retention fractions for single and binary populations, 
all NSs. For comparison, we show the total retention fraction of a binary population with the 
Arzoumanian et al. (2002) kick distribution (long-dashed line). 

2. Production and retention 

In our studies we adopt that a NS can be formed as a result of either a core-collapse 
(CC) supernova or an electron capture supernova (ECS). In the latter case, three possi- 
bilities are considered: ECS during normal stellar evolution, accretion induced collapse 
(AIC) and merger induced collapse. For the case of CC NSs, we adopt that a supernova 
was accompanied by a natal kick in accordance with the distribution by Hobbs et al. 
(2005). For ECS NSs, we adopt that the accompanying natal kick is 10 times smaller. 
We find that even considering a stellar population with 100% primordial binaries, the 
retention fraction of CC NSs is very small (Fig. 1) and the resulting number of retained 
NSs is just a few per typical dense globular cluster of 2 x 10''' Mq. NSs formed via dif- 
ferent ECSs channels are retained in reasonable numbers, providing about 200 retained 
NSs per typical GC, or more than a thousand in a cluster like 47 Tuc (similar numbers 
were found also in Kuranov & Postnov, 2006). Therefore, in contrast to the population 
of NSs in the Galaxy, the population of NSs in GCs is mainly low-mass NSs made by 
ECS. 



3. X-ray binaries 

In our simulations, we find that a typical GC can contain up to 2 LMXBs with a MS 
companion (most likely observed, at any particular time, as qLMXBs) and up to one 
LMXB with a WD companion (ultra-compact X-ray binaries, UCXBs). The scatter in 
the average number of observed LMXBs per cluster in independent simulations is very 
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Figure 2. The collision numbers 7 (Verbunt & Hut 1987) and numbers of LMXBs in simulated 
clusters. The solid line corresponds to nLMXB ~ (2.7±6) +O.O287, 7 is per W^Mq, as in Pooley 
& Hut (2006). The error bars correspond to the scatter in our simulations. The standard model 
has a central core density of 10® pc~"^, velocity dispersion 10 km/s and half-mass relaxation time 

I Gyr. The low-p, med-p, standard and high-p models have central core densities 10'^, 10*, 10® 
and 10® pc""', respectively. Our low-a model has velocity dispersion 5 km/s, while long-irh has 
half-mass relaxation time 3 Gyr. 

large - e.g., for UCXBs, it can vary between 0.1 and 1.1. In the case of Terzan 5 and 47 
Tucanae, the average number of LMXBs formed per Gyr, at the age of 11 Gyr, is '^ 5 for 
NS-MS LMXBs and ~ 8 for UCXBs. These numbers are in general agreement with the 
observations. Overall the numbers of NS that gain mass via mass transfer (MT) through 

II Gyr of cluster evolution are high: for our 47 Tuc model, about 40 NS-MS binaries and 
more than 70 UCXBs. As we observe fewer MSPs in these GCs, while the rate of LMXB 
formation in simulations is consistent with the observations, we conclude that not all 
NSs that gain mass via MT become currently active MSPs. 

We analyzed how the specific number of LMXBs hlmxh in our simulations depends 
on the specific collision frequency (see Fig. 2). For the case when only core density is 
varied, hlmxe depends linearly on 7. Variation of other cluster dynamical properties 
leads to deviation from such a linear dependence, which may explain the scatter in 7 in 
the observed GCs. 



4. Millisecond pulsars 

Suppose that all mass-gaining events in the life of a NS - mass transfer, physical 
collision with a red giant, common envelope hyper- accretion or merger - can lead to NS 
recycling. In this case we found that as many as 250 and 320 potential MSPs are made 
in our simulations of clusters like Terzan 5 or 47 Tuc, accordingly (the corresponding 
numbers of retained NSs are ~ 500 and ^ 1100). Although these numbers correlate 
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Figure 3. bMSPs in simulated models of 47 Tuc and Terzan 5 compared to observed bMSPs. The 
simulated populations correspond to several independent runs and represent a larger population 
than in the observed clusters. Observed bMSPs are shown with circles; triangles - bMSP formed 
via binary exchanges; stars - via tidal captures; squares - via physical collisions; diamond - 
primordial binaries. Cross signs mark eccentric bMSPs (e ^ 0.05) and solid symbols mark 
systems with a non-degenerate companion (in the case of simulations) or observed eclipsing 
systems. 



well with the formation rate of LMXBs, it greatly exceeds the numbers of observed and 
inferred MSPs in both clusters, which are 33 (in Terzan 5; perhaps 60 total) and 22 (in 
47 Tucanae; perhaps 30 total). 

We analyzed the population of NSs that gained mass. We found that bMSPs formed 
from primordial binaries, where a common envelope event led to AIC, create a population 
of potential bMSPs with relatively heavy companions, in circular orbits with periods from 
one day to several hundreds of days. This population is not seen in either Terzan 5 or 47 
Tuc. We considered primordial binaries that evolved through mass transfer from a giant 
donor after a NS was formed via AIC. Even though bMSPs that have similar periods, 
companion masses and eccentricities are present in Terzan 5, there are no such systems 
in 47 Tucanae. Also, bMSPs made from primordial binaries after AIC must inevitably be 
formed in low-dense clusters, but no bMSPs are observed there. These facts tell us that 
either AIC does not work, or the kicks in the case of AIC are stronger then we adopted, 
or a NS formed via AIC has such a strong magnetic field, that surface accretion does not 
occur. 

A common understanding of MSP formation is that the NS is recycled through disk 
accretion, where a NS is spun up only if the accretion rate is not too low, Af ^ 3 x 
lO^^AfEdd, where AfEdd is the Eddington limit (for a review see, e.g.. Lamb & Yu 2005). 
In a UCXB, soon after the start of mass transfer, the accretion rate drops very quickly. 
After 1 Gyr, it is less than lO^^AfEdd- Such a MT leads to a spin-down of the previously 
spun-up NS, and no MSP is formed. Support for this statement is given by the fact that 
no UCXBs (those that have WD companions) are visible as MSPs (Lamb & Yu 2005). 

The requirement of steady spin-up through disk accretion implies that not all physical 
collisions will lead necessarily to NS spin-up. In the case of a physical collision with a 
giant, the NS will retain a fraction of the giant envelope, with a mass of a few hundredths 
of A^o (Lombard! et al. 2006). Immediately after the collision, this material has angular 
momentum and most likely will form a disk. We adopted therefore that in the case of 
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a physical collision with a giant, an MSP can be formed, but, in the case of any other 
physical collision, the NS will not be recycled. 

Considering all the exclusions described above, we form in our simulations at least 
15 ± 7 MSPs for Terzan 5 and 25 ± 4 MSPs for 47 Tuc (for the formed population of 
bMSPs, see Fig. 3). The values for Terzan 5 are somewhat uncertain due to uncertainty 
in the properties of this heavily reddened cluster. The total number of NSs that gain 
mass in the simulations by one or other way are 250 and 320 in Terzan 5 or 47 Tuc, 
accordingly. 



5. Conclusions 

We studied the formation and evolution of NSs in globular clusters. We find that NS 
formation via different channels of ECS is very effective in GCs and provides most of the 
retained NSs. Having as many as a few hundred retained NSs per typical GC, or about 
1100 per cluster like 47 Tuc, produces LMXBs in numbers comparable to observations. 
We note that if AIC does not lead to the formation of NSs, then the number of formed 
NSs is reduced only by ~ 20%, but the number of appearing LMXBs is decreased by 2-3 
times (per Gyr, at the cluster age of 11 Gyr), although it may still be consistent with the 
observations, given the large scatter in the simulations. We find that up to half of NSs 
could gain mass after their formation through mass transfer, hyper-accretion during a 
common envelope, or physical collision. It is likely that most of these mass-gaining events 
do not lead to NS spin-up, and that only a few per cent of all NSs appear eventually as 
MSPs, implying that there is a large underlying population of unseen NSs in GCs. 
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